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Abstract 

Torrefaction is a thermal treatment step in the relatively low temperature range of 225-300 °C, which aims to produce a fuel with increased 
energy density by decomposing the reactive hemicellulose fraction. The weight loss kinetics for torrefaction of willow, a deciduous wood type, was 
studied by isothermal thermogravimetry. A two-step reaction in series model was found to give an accurate description. For the two steps, 
activation energies of 76.0 and 151.7 kJ/mol, respectively, and pre-exponential factors of 2.48 x 10 4 and 1.10 x 10 10 kg kg -1 s -1 , respectively, 
were found. The first reaction step has a high solid yield (70-88 wt%, decreasing with temperature), whereas less mass is conserved in the second 
step (41 wt%). The fast initial step may be representative for hemicellulose decomposition, whereas the slower subsequent reaction represents 
cellulose decomposition and secondary charring of hemicellulose fragments. The kinetic model is applied to give recommendations for industrial 
torrefaction process conditions, notably operating temperature, residence time and particle size. 

© 2006 Elsevier B.V. All rights reserved. 

Keywords: Thermo gravimetry; Torrefaction; Pyrolysis; Kinetics; Wood; Hemicellulose; Cellulose; Lignocellulosic fuels 


1. Introduction 

Thermal degradation of wood is a complex research topic 
because wood contains different fractions. On microscopic 
scale, wood cells are composed of so-called microfibrils, 
bundles of cellulose molecules ‘coated’ with hemicellulose. In 
between the microfibrils (and sometimes within the amorphous 
regions of the microfibril) lignin is deposited. These wood 
fractions show different thermal behaviour. Three zones may 
therefore be distinguished in weight loss curves of wood: 
hemicelluloses, the most reactive compounds, decompose at 
temperature in the range of 225-325 °C, cellulose at 305- 
375 °C and lignin gradually over the temperature range of 250- 
500 °C [1]. This research addresses specifically the first zone, 
i.e. decomposition of hemicellulose in the relatively low 
temperature range of 225-300 °C, a process also known as 
torrefaction [2-5]. The solid product of such a process, roasted 
or so-called torrefied wood, has found applications as a 
barbecue fuel and firelighter [4] . It may be attractive in future as 
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a renewable fuel with increased energy density (compared to 
raw wood) for gasification and/or co-combustion. 

Although a wealth of research data is available about pyrolysis 
of wood, relatively few address the temperature range of 
torrefaction. Therefore, more research is required to obtain better 
understanding of the chemical and physical processes occurring 
in torrefaction, especially reaction mechanism, reaction kinetics, 
optimization of process conditions and selection of particle size. 
This article focuses on the determination of the weight loss 
kinetics of wood torrefaction by experiments in a thermogravi- 
metric analyser (TGA). In a second article, the products formed 
in the torrefaction process are quantified and analysed, with 
emphasis on the properties of torrefied wood. 

The weight loss kinetics of lignocellulosic fuels have been 
determined in many different experimental devices. Such 
devices include thermogravimetric analysers, tube furnaces and 
fluidized bed reactors. Experiments can be carried out under 
dynamic conditions, in which the sample is submitted to an 
assigned heating rate, but experimental heating rates are often 
much slower than those achieved in process equipment such as 
gasifiers or combustors. The research described below uses 
static conditions, in which the sample is kept at a constant 
temperature. In this case, the heating rate is important too, 
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because the sample has to be heated from environmental 
temperature to the reaction temperature at which the weight 
loss kinetics is studied. If the heating rate is too high, the results 
may be affected by heat transfer limitations within the sample. 
If a too slow heating rate is applied, the weight loss that takes 
places during the warm-up phase is not negligible, which 
complicates the interpretation and deduction of kinetic data. 
However, for the low temperatures applied in torrefaction (225- 
300 °C), the warm-up phase is relatively short, even for low 
heating rates of 10-100 K/min common in TGA equipment and 
the error introduced by the weight loss during this phase is 
acceptable. 

An important research parameter is the type of biomass, the 
composition of which determines its behaviour in the torrefaction 
process. Deciduous wood types as well as coniferous wood types 
are considered. These types are often referred to as hard wood 
and soft wood, respectively. There is not much difference in the 
distribution of the three wood fractions between these types. 
Coniferous wood may contain slightly more lignin than 
deciduous wood (25-35 wt% versus 18-25 wt%), slightly less 
cellulose (35-50 wt% versus 40-50 wt%) and, on average, 
comparable amounts of the hemicelluloses (20-32 wt% versus 
15-35 wt%) [6] . However, the composition of the polysugars that 
form the hemicelluloses fraction is very different. In the case of 
deciduous wood, the hemicelluloses contain 80-90 wt% of 4-O- 
methyl glucoronoxylan (referred to as xylan) whereas they 
contain 60-70 wt% of glucomannan and 15-30 wt% of 
arabinogalactan for coniferous wood [6] . The thermal behaviour 
of these components may be different. Therefore, it is an 
important question whether the thermal behaviour of deciduous 
and coniferous wood types is the same, i.e. do they react with 
comparable rates in the torrefaction temperature range. 
Furthermore, weight loss kinetics of wood is compared to wood 
components such as xylan. 

Finally, the kinetic model obtained in this work is applied to 
make recommendations for industrial torrefaction process 
conditions, such as temperature, residence time and particle size. 

2. Experimental 

2.7. Materials 

Several biomass types have been tested in the experiments: 
deciduous wood (beech and willow), coniferous wood (larch) 
and straw. The deciduous wood types were obtained from 
Rettenmaier, Germany, the coniferous wood from Praxis, 
Netherlands and wheat straw was taken from a field in North- 
Holland. Particle sizes were in the range of 0.7-2.0 mm in all 
cases, except for straw where it was <5 mm. 

Model compounds of wood, such as cellulose and 4-O-methyl 
glucuronoxylan, extracted from oak spelt, were also used. These 
were purchased from Sigma-Aldrich in powdered form. 

2.2. Equipment and procedures 

TGA measurements were carried out with a Perkin-Elmer 
Pyris 6 TG apparatus with autosampler. The analysed sample 


weights varied from 2 to 10 mg. Analysis was carried out with 
a nitrogen flow rate of 20 ml (STP)/min. Ceramic crucibles 
loaded with biomass particles were placed inside the TG 
apparatus where the weight was constantly measured. The 
temperature was varied from 230 to 300 °C and the heating 
rate from 10 to 100 °C min -1 . A temperature program 
consisting of a dynamic heating period followed by an 
isothermal heating period was applied. During the dynamic 
heating period the biomass was thoroughly dried. The weight 
loss graphs of the dry biomass were used for determination of 
reaction kinetics. 

Due to variations in ash content of the various samples, the 
measured weights were corrected in the following way: 

(K\ = WtgA ~ ^ash 

UoAxp initial - W ash () 

where TVi nitia i is the solid weight of the dry sample, W ash the 
weight of the ash in the dry sample and W TGA is the solid weight 
that is measured as a function of time. 

3. Results and discussion 

3.1. Thermal degradation of different biomass types 

Fig. 1 shows the weight loss curves of the various biomass 
types at 248 and 267 °C obtained from isothermal TGA 
experiments. The weight loss observed during heating of the 
sample from 200 °C, the temperature at which thermal 
decomposition begins to occur, to the required temperature 
is relatively small, except for xylan. From Fig. 1 it can be 
concluded that xylan, the main hemicellulose component of 
deciduous wood, is the most reactive component. It starts 
decomposing around 200 °C and has a high weight loss. The 
cellulose decomposition rate is very low in the temperature 
range used. This is in agreement with results of other studies, 
e.g. Gaur and Reed [7] mention that cellulose may 
depolymerise at low temperatures, but volatilization of the 
cellulose is very slow below 250 °C. At 267 °C, limited weight 
loss of cellulose is found. 

The hardwoods beech and willow have comparable 
reactivity. The weight loss observed lies between that of 
xylan and cellulose, which was expected as wood contains 
these fractions (and other fractions, such as lignin). High xylan 
content also explains the relatively high weight loss of wheat 
straw, although catalytic effects due to the presence of mineral 
matter could also play a role. Finally, the coniferous larch 
reacts a lot slower than the hardwoods. Although some 
researchers have attributed the lower weight loss for 
coniferous wood to higher lignin content, Ramiah [8] has 
experimentally shown that differences in composition of the 
hemicellulose fractions in deciduous and coniferous wood 
(xylan- and mannan-based, respectively) may be the main 
explanation. 

As a preliminary conclusion, hardwood looses considerably 
more weight than softwood and therefore a higher increase in 
energy density (J/g) may be expected, a topic that is studied in 
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( a ) Time (s) 



(b) Time (s) 

Fig. 1. TGA of various biomass compounds, (a) at 248 °C and (b) at 267 °C. 
Heating rate 10 °C min -1 , particle size 0.5-2 mm; dotted line is the heating 
curve. 


with two successive reactions, as given in Eq. (2). In such a 
model, xylan is assumed to form an intermediate reaction 
product, which is a solid with a reduced degree of 
polymerization. This intermediate then reacts to the final 
product. The first reaction was found to be substantially faster 
than the second one. 


>)>iB + (1 -yi)Vi 

(2a) 

•yiC + (1 — yi)V 2 

(2b) 


The yields of solid product are denoted as the parameters y i and 
y 2 for the first and second reactions, respectively. Varhegyi et al. 
have regarded these parameters as temperature-independent, 
e.g. for xylan decomposition, constant values for the solid 
yields of the first and second reactions were reported. There¬ 
fore, the model fails to take into account the change of final char 
yield (which equals the product of yi and y 2 ) with temperature 
[13]. Nevertheless, the final char yield can be expected to 
decrease with temperature, as the hemicellulose component 
is dehydrated to a larger extent. 

Di Blasi and Lanzetta [13] used isothermal TGA measure¬ 
ments for evaluation of the kinetic parameters of torrefaction of 
xylan. They modeled the kinetics by a combination of a two- 
step mechanism with parallel reactions for the formation of 
solids and volatiles. 



A-»B-*> C 


the second paper. This makes hardwood an attractive feedstock 
for torrefaction processes. In the rest of this section, a kinetic 
model is therefore developed for torrefaction of willow, a 
classical hardwood. 


The solid yields are given by: 

_ h 
k\ + k Vl 


3.2. Kinetic model for willow torrefaction 


yi 


k2 

k 2 + ky 2 


(4a) 

(4b) 


Because the pyrolysis of wood is a very complex set of 
reactions, simplified kinetic models are generally applied. For 
pyrolysis of cellulose and biomass at temperatures higher than 
300 °C, the Broido-Shafizadeh model [9,10] is the mostly used, 
consisting of an activation step followed by two subsequent 
reactions. Other models consisting of three first-order 
consecutive reactions or of a simple first-order reaction (global 
one step) have been proposed as well [11,12]. For hemi¬ 
cellulose decomposition, a global one-step reaction gives a very 
rough impression of kinetics of the overall reaction but it is 
difficult to compare parameters consistently with other 
researchers [13]. Orfao et al. [14] concluded that hemicellulose 
decomposition could not be modeled by simple kinetics. 
Research done on multi-step reaction mechanisms gives better 
results [11,12,15,16]. Varhegyi et al. [11,12] have studied the 
reaction kinetics of thermal decomposition of xylan by 
dynamic TGA-MS. They found a model of competitive 
reactions could not describe their data, and proposed a model 


When the activation energy of k V] (respectively, ky 2 ) is higher 
than that of ki (respectively, k 2 ), the solid yield drops with 
increasing temperature. So the introduction of separate forma¬ 
tion of different product classes, which may be questioned from 
the viewpoint of analytical chemistry [17], is a mathematical 
approximation to describe the fact that the ratio between 
formed solids and volatiles depends on temperature. 

For this research, it was decided to verify whether a two-step 
mechanism (as given in Eq. (3)) is also suitable to describe the 
weight loss kinetics of wood in the torrefaction temperature 
range. This study does consider solid formation and volatile 
formation as parallel reactions and aims to determine k\, k 2 , ky x 
and ky 2 . Our aim was to develop a simple model that describes 
the weight loss of willow wood in the temperature range of 
torrefaction—without exact knowledge of the wood composi¬ 
tion. Other wood types or biomass types with different 
composition will have different char yields, and also different 
‘kinetic’ parameters. Therefore, the char yields and reaction 
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rates of the first and second step for willow wood are compared 
with literature data for xylan. 

The differential rate equations are given for the solids by 
Eqs. (5a)-(5c). 

rA = = ~( ki + ( 5a ) 

re = = kx [A]" A - (k 2 + *y 2 )[Bp (5b) 


rc = d ^ = k 2 [Br (5c) 

Similar equations may be derived for the volatile products. 

If all reactions are assumed to be of first order, the system of 
equations can be solved analytically. Integration of the 
differential equations, with the initial condition that only A 
is present at the beginning of the reactions, and addition gives 
an expression for the relative solid weight: 

M=[ A] + [B] + [C] (6) 


M t _ / I" kjKi - k x k 2 
M 0 ~\ [ K\ (K2 K \) 

k\k2 

where 


e 


-Kit 


+ 


— k\K2 + k\k2 
K 2 (K 2 -K l ) 


e 


-K 2 t 


(7) 


K\ — k\ -\- ky ] (8a) 

K 2 = k 2 ky 2 (8b) 


The final char yield, which can be determined experimen¬ 
tally, is equivalent to the last term in Eq. (7): 


Ct —> oo k\k2 

M 0 _ K X K 2 


(9) 



Fig. 2. Experimental weight loss curve for willow decomposition at various 
final temperatures. Starting weight defined at 200 °C; time includes a heating 
period of 3-10 min (heating rate = 10 °C min -1 ). 


reaction, as is the case for xylan decomposition, a demarcation 
time exists after which the first reaction is completed and the 
second one starts. For decomposition of willow wood, it was 
possible to establish a demarcation time at temperatures in 
excess of 250 °C. From this temperature, it seems that 
secondary charring reactions begin to occur. 

To determine the reaction orders, a plot of the logarithm of 
the non-dimensional solid weight versus time is shown in Fig. 3 
(for clarity, experimental data points are shown with 5-10 min 
intervals). Straight lines are observed before and after the 
demarcation time, so that a model with two first-order reactions 
is valid. This is advantageous as it gives an analytical 
expression for the weight loss, whereas heat transfer limitation 
calculations become a lot more complicated for higher order 
models. 

The kinetic parameters could be determined by a graphical 
method or by a numerical method. In the graphical method, a 
demarcation time t should be determined at each experimental 
temperature. The sample weight at this time is assumed to 
consist only of intermediate pseudo-component B (as A has 
reacted away and C has not formed yet), i.e. B . However, for 


For higher reaction orders, e.g. a second-order reaction fol¬ 
lowed by a first-order sequential reaction, the system of equa¬ 
tions is non-linear and of stiff nature. This system can be solved 
numerically using a suitable algorithm, e.g. in Matlab Release 
12, a one-step solver based on a modified Rosenbrock formula 
of order two may be applied. 

3.3. Fitting of weight loss curves 

Fig. 2 shows the experimental weight loss curve for 
decomposing willow at various temperatures for a time interval 
of 3 h. This interval includes a heating period of 3-10 min, 
required to heat the sample from 200 °C to the final temperature 
(at a heating rate of 10 °C min -1 ). The weight loss curves for 
wood and wood components are modeled by a reaction in series 
model as this gives more accurate results than a global one-step 
reaction. If the first reaction is must faster than the second 


Time [min] 



Fig. 3. Logarithmic relative weight of willow vs. time, at various final 
temperatures. Starting weight defined at 200 °C; time includes a heating period 
of 3-10 min (heating rate = 10 °C min -1 ). 
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our data, Fig. 3 shows that an exact demarcation time is difficult 
to establish, due to overlapping of the reaction stages. 
Therefore, a numerical approach is more accurate than a 
graphical approach, and the best fit for all kinetic parameters 
was achieved in Matlab by minimizing the sum of squares 
function: 



3.4. Kinetic parameters 


The kinetic parameters obtained using the method described 
are as follows: 


k\ = 2.48 x 10 4 exp 


kvt = 3.23 x 10 7 exp 
k 2 = 1.10 x 10 10 exp 


ky 0 — 1.45&2 


-75976 \ 

(Ha) 

RT J 

— 114214 \ 

(Hb) 

RT J 

— 151711 \ 

(11c) 

RT ) 


(Hd) 


All rate constants are expressed in kg kg -1 s -1 , the temperature 
T in K and the gas constant R in J mol -1 K -1 . The kinetic 
parameters have been obtained using experiments with a low 
heating rate of 10 K/min, taking into account the weight loss 
during the heating phase. 

In agreement with literature, it was found that the first 
reaction stage is remarkably faster than the second stage. 
Activation energies for the first and second stages are very close 
to those reported by Branca and Di Blasi [18] of 76 and 143 kJ/ 
mol, respectively. The latter values have been determined by 
isothermal thermogravimetry for beech wood, for which the 
thermal behaviour has been demonstrated to be very similar to 
willow wood. Compared to the parameters published by Di 
Blasi and Lanzetta for xylan torrefaction, willow wood 
decomposes approximately one order of magnitude slower 
than xylan. This can be partly explained because willow 
contains only 25-30 wt% of xylan, but also because extraction 
of pure xylan from a wood matrix may make the xylan more 
reactive. 


3.5. Model verification 

To verify the final char yield predicted by the model, TGA 
experiments with a longer reaction time of 435 min were 
carried out (including a heating period of 7-10 min at a 
heating rate of 10 °C min -1 ). The results are shown in Fig. 4. 
At temperatures of 280-300 °C, predicted char yields of 29- 
31% showed good agreement with experimental data, 
whereas determination of the final char yield at 270 °C 
was not possible because the weight had not stabilized yet 
within 435 min. Careful consideration of the experimental 
weight curves teaches that these do not become completely 



Fig. 4. Experimental and modeled relative weight of willow vs. time, for 
various temperatures. Starting weight defined at 200 °C; time includes a heating 
period of 7-10 min (heating rate = 10 °C min -1 ). 


horizontal, but the weight continues to decrease very slowly. 
Even when the reactive components (hemicellulose and 
cellulose) have reacted away, the remaining lignin-rich 
fragment continues to react. To take this phenomenon into 
account, the model of Eq. (3) should be extended, which 
could be done in two possible ways. One option would be to 
split pseudo-component A into two fractions: a lignin-rich 
fraction and a fraction of labile components, so that the 
reaction in series model is extended with a parallel reaction 
describing lignin decomposition. Alternatively, a third 
reaction may be added to the reaction in series model, as 
in Branca and Di Blasi [18]. For the low temperature range 
relevant in this paper, however, a two reaction in series model 
is sufficiently accurate. 

The final char yield is the product of the solid yield for the 
first decomposition reaction (>']) and the second decomposition 
reaction (y 2 )- For the temperature range studied, y x decreases 
from 88% at 230 °C to 70% at 300 °C, whereas y 2 is fairly 
constant at 41%. The first reaction with a relatively low weight 
loss could very well be representative of decomposition of the 
reactive xylan in willow wood, since hardwood contains up to 
30% of xylan. The high weight loss in the second reaction can 
be explained by decomposition of the other fractions contained 
in the wood, notably the cellulose fraction, perhaps in 
combination with charring of the remaining hemicellulose 
fraction. This is also plausible because the activation energy 
found for the second weight loss step is close to literature values 
[19]. It may be that the primary decomposition products of 
xylan decomposition, e.g. acids, initiate decomposition of 
cellulose fraction. This hypothesis has been formulated 
previously [5] and contradicts the approach of others who 
have considered the rate of biomass decomposition as the sum 
of the rates of main biomass components that react individually 
[14,20]. 

The model has also been tested at higher heating rates up to 
100 °C min -1 . Fig. 5 shows the sample temperature and relative 
weight as a function of reaction time for heating rates of 10, 50 
and 100 °C min -1 and a final temperature of 260 °C. It can be 
observed that the trend in weight loss is predicted accurately 
also at higher heating rates. At a heating rate of 100 °C min -1 , 
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Time [min] 


-10 °C/min 

* exp (10 °C/min) 

-50 ®C/min 

o exp (50 °C/min) 

.100°C/min 

+ exp (100 °C/min) 

-temperature (10 °C/min) 

-temperature (50 °C/min) 

.temperature (100 °C/min) 


Fig. 5. Experimental and modeled relative weight of willow vs. time, with varying heating rates and final temperature of 260 °C. Bottom lines indicate sample 
temperatures. 


the required heating time is less than a minute and the 
associated weight loss less than 1%. 

3.6. Recommendations for industrial torrefaction process 
conditions 

The kinetic data may be used to make recommendations for 
industrial torrefaction process conditions and for process scale- 
up. In torrefaction, the fast initial reaction, in which the highly 
reactive hemicellulose is decomposed, is important. This 
reaction stage should go more or less to completion. The second 
reaction stage, which could represent cellulose decomposition 
and secondary charring reactions of hemicellulose reaction 
products, takes much more time to complete, which may lead to 
uneconomically large equipment. Fig. 6 shows the time 
required to obtain the maximum yield of intermediate product 
B. This time reduces from over 3 h at 230 °C to 10 min at 
300 °C. The latter reaction time appears short enough for 
commercially viable processes, since the estimated reactor 
volume may be fitted into a single reactor even at large 
capacities up to 10 kg s _1 . 



Fig. 6. Time required for maximum yield of intermediate product B and 
corresponding yields of A, B and C at different temperatures. 


The maximum yield of B, which is also shown in Fig. 6, 
decreases with temperature since more volatiles are formed at 
higher temperature. This may be favourable if these volatiles 
consist of species that carry little energy, such as water vapour 
or carbon dioxide. However, at temperatures above 300- 
320 °C, fast thermal cracking of cellulose may cause tar 
formation, so that operation below 300 °C is recommended. 
This also ensures a good selectivity to the intermediate product; 
the selectivity to B decreases with temperature because the 
activation energy of the secondary reactions is higher than that 
of the primary reactions. 

Apart from temperature and residence time, another 
important process parameter for torrefaction is the particle 
size. To ensure the minimum residence time, it is important that 
the rate of heat transfer to and within the particle is fast, as 
compared with the reaction rate. This means that the solid 
temperature will be essentially the same as the surrounding 
temperature, and the overall controlling factor is the intrinsic 
kinetics. The ratio between heat convection and conduction 
rates is expressed by the Biot number: 

Bi = (12a) 

A 


while the ratio of heat convection rate and reaction rate is 
expressed by the (external) Pyrolysis number [21]: 


Py' = 


a 

K x pc v r v 


(12b) 


In these equations, a is the external heat transfer coefficient in 
W m -2 K -1 , r p the radius of the particle in m (assuming 
spherical particles), A the thermal conductivity in W m 1 K -1 , 
p the density in kg m “ 3 and c p is the heat capacity in J kg” 1 K -1 of 
the biomass particle. 

The desired kinetic control is achieved when the following 
criteria for absence of heat transfer limitations are met: (a) Biot 
number much smaller than unity, so that heat conduction within 
the particle is much faster than heat convection to the particle 
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Biot number [-] 


Fig. 7. Biot number and Pyrolysis number (for reaction at 250 and 300 °C) as a 
function of biomass particle diameter. Biomass properties: 1 = 0.1 W m -1 K _1 , 
p = 1500 kg m -3 and c p = 2000 J kg -1 K -1 . External heat transfer coefficient 
a = 100 W m -2 K“\ 

and (b) external Pyrolysis number much larger than unity, 
which indicates that heat convection to the particle is much 
faster than the chemical reactions taking place. Fig. 7 shows 
these parameters at reaction temperatures of 250 and 300 °C for 
an external heat transfer coefficient of 100 W m -2 K 1 . At a 
Biot number of unity, equivalent to a particle size of 2 mm, the 
Pyrolysis number is sufficiently high. For pyrolysis of wood 
below 300 °C, it follows that the reactions are the rate limiting 
step, provided that the particle size is sufficiently small and 
internal heat transfer does not play a role. 

4. Conclusions 

The kinetics of torrefaction reactions in the temperature 
range of 230-300 °C can be described accurately by a two-step 
mechanism. The first step is much faster than the second step, 
so that these steps can be demarcated in time. The first step is 
representative of hemicellulose decomposition, while the 
second step represents cellulose decomposition. The solid 
yield for the first step is higher than for the second step: 70-88% 
(decreasing with temperature) versus 41%. This may be 
explained because deciduous wood, such as willow, contains 
less xylan (the reactive component in its hemicellulose fraction) 
than cellulose. 

Comparison of the kinetic parameters of willow with those of 
xylan shows that xylan reacts approximately one order of 
magnitude faster than willow, a classical hardwood containing up 
to 30% xylan. It is recommended to determine reaction kinetics 
also for softwood, such as larch. However, initial experiments 
have shown that softwood is much less reactive and therefore 
lower rate constants can be expected. Research into thermal 


degradation of model components such as glucomannan, which 
constitutes 60-70 wt% of the hemicellulose fraction in soft¬ 
woods, is recommended to confirm this finding. 
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